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Abstract 

Motivated by the description of the graphene electronic structure in terms of the relativistic Dirac equa- 
tion, a generalization to four dimensions yields a strictly local fermion action describing two species and 
possessing an exact chiral symmetry. This is the minimum number of species required by the well known 
"no-go" theorems. 
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The structure of graphene, a single layer of graphite consisting of a hexagonal lattice of carbon 
atoms, has attracted considerable attention recently both from the experimental front and the fact 
that the low electronic excitations are described by the Dirac equation for massless fermions yj]. 
From the point of view of a particle physicist, this structure has two particularly striking features. 
First, the massless structures are robust for topological reasons related to chiral symmetry. Second, 
it achieves this symmetry in a manner that involves the minimum number of effective massless 



fermions required by the famous "no-go theorems" for lattice chiral symmetry Q3|, |4fl. 

Given the importance of chiral symmetry in particle physics and the difficulties with imple- 
menting it with a lattice regularization [5], it is natural to ask whether these properties of the 
graphene electronic structure can be extended to four dimensions. Indeed, this is possible, and 
provides a remarkable fermion action with an exact chiral symmetry and manifesting two species 
of massless states, the minimal number consistent with chiral symmetry. This action is strictly 
local and thus will be vastly faster in simulations than either the overlap operator or domain wall 
fermions, the only other known ways to have chiral symmetry with only two flavors. 

Recently a chiral gauge theory structure on two dimensional graphene has been proposed \^\. 
Given that we do not yet have a lattice regularization of the standard model, it would be particularly 
interesting if this construction could be extended to the four dimensional lattices presented here. 

Although well known, it is useful to review briefly the standard two-dimensional graphene 
band structure. We will closely parallel this derivation for the four-dimensional case. Ignored 
here are all but the pi orbitals in a tight binding approximation. We also ignore the electron 
spin, each component of which gives an equivalent structure. Our electrons hop from neighbor to 
neighbor around a fixed underlying hexagonal lattice. A fortuitous choice of coordinates makes 
the problem straightforward to solve. As sketched in Fig. (OQ), orient a graphene surface with one 
third of the bonds horizontal, one third sloping up at 60 degrees, and one third sloping down. It 
is then convenient to collapse the atoms at the opposite ends of each horizontal bond together and 
call this unit a lattice "site," as enclosed in ellipses in the figure. For each site, let denote the 
creation operator for an electron on the left atom, and correspondingly let create an electron on 
the right atom. The commutation relations are the usual 

[a XuX2 ,al^] + = [b Xl ^,b^] + = 8 XlA 8 X2A (1) 

with the a type operators anti-commuting with the b's. Finally, it is useful to label the sites using 
a non-orthogonal coordinate system with axes x\ sloping up at 30 degrees intersecting the corre- 



FIG. 1: Organize the graphene structure into two-atom "sites" involving horizontal bonds as shown by 
ellipses. Call the left hand atom of each site type a and the right hand atom type b. The coordinates of the 
sites are labeled along the non-orthogonal x\ and a"2 axes. 



sponding sites, and similarly X2 sloping down at 30 degrees. All of this is illustrated in Fig. ©. 

With these conventions, the Hamiltonian of interest involves only nearest neighbor hoppings 
between a and b type sites 

H = ,x 2 ( a xi ,x 2 b*i ,*2 + bx\ ,x 2 a xi ,*2 

+ a l 1 + l,x 2 b ^2+ b l 1 -l,x 2 a ^2 (2) 

+ a !i ,x 2 - 1 ,x 2 + b Xi X2+l a Xl tX2 ) . 
Here K is the basic hopping parameter. The phase of K is a convention; here I consider positive 
real hopping. To diagonalize this Hamiltonian go to momentum space 



1, 2 L„ 2% 2% PUP2 



This brings the Hamiltonian to the simple form 

H = K S-&%t al upi b PuP2 {l+e-^+e^) 
+bl l , P2 ZpuP2(l+e +m +e- ip2 : 
The problem reduces to diagonalizing a two by two matrix of form 



(4) 



H{p uP2 )=K \ ° ■ I (5) 




with 

z=l + e~ ipi +e +ip2 . (6) 
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The energy eigenvalues are 



E{pupi) = ±K\z 



(7) 



From Eq. © it is easy to see that the energy vanishes only at two points, p\ = P2 = ±2n/3. 

These are known as "Fermi points" (2D and their robustness can be seen by considering contours 
of constant energy. These are closed curves of constant \z\ in pi,p2 space . The important point 
is that for such a contour near one of the zero energy solutions, the phase of z wraps non-trivially 
around a circle. This non-contractable mapping indicates that on reducing the energy and shrinking 
the curve to a point, the magnitude of the energy at this point must vanish. This is the mechanism 
that prevents a band gap from opening in the spectrum. 

This robustness is associated with a chiral symmetry. Because the hoppings are always between 
a and b type sites, we can change the sign of the energy by taking b — > —b. This is equivalent to the 
statement that G3 anti-commutes with the Hamiltonian. For the four-dimensional generalization, 
this will become the anti-commutation of 75 with the Dirac operator. 

We wish to extend this formalism to the four-dimensional case. We want an operator D to 
insert into the Euclidian path integral via the fermion action \j7D\]/. For low energy excitations this 
operator should reduce to two massless Dirac fermions, and this reduction should be robust due to 
a chiral symmetry. Essentially all Dirac operators used in practice for lattice gauge theory satisfy 
what is called a "75 Hermiticity" condition 



where 75 is the usual four by Dirac matrix. It is Hermitean and squares to the unit matrix. A 
specific realization of this and the other Dirac matrices will appear below. Using this we can 
construct a Hermitean "Hamiltonian" 



with which we will parallel the two dimensional discussion. It is important to remember that this 
is not the Hamiltonian of the three dimensional quantum system, but a convenient operator for 
leading us back to D. In four-dimensional space, the analog of the curves of constant "energy" 
are three dimensional manifolds. To maintain a topological argument in analogy to the two di- 
mensional case, we want to consider the situation where these surfaces wrap non-trivially around 
a three sphere, an S3. For this purpose it is quite natural to maintain the form of Eq. ©, but extend 
Z to two by two matrices in a quaternionic space. That is, take 



y 5 Dy 5 =D 



(8) 



H = y 5 D 



(9) 



z = ao + ia- o 



(10) 



4 



with a ^ a real four vector and o denotes the traditional Pauli matrices. With z being a two by 
two matrix inserted into the two by two matrix of Eq. ©, we wind up with a four by four matrix, 
the same dimension as used in the usual Dirac equation. As before, vanishing energy states occur 
when z vanishes, which now corresponds to vanishing as a four vector. The goal is to construct 
our Hamiltonian so that that constant energy surfaces that wrap around zero energy points within 
the Brillouin zone will involve a non-trivial mapping in the quaternionic space. Because of the 
periodicity of the Brillouin zone, these zero energy points must appear in pairs so that the overall 
wrapping will vanish. Indeed, this is the famous no-go theorem [[3i|4|]. 

We want a construction giving precisely one and only one such pair. We also want to involve 
only local couplings, i.e. with only simple trigonometric functions of the momenta appearing 
in the dispersion relations. Because of the robustness of the zeros, if such a construction exists, 
it is clearly not unique. To find one such solution, start with a regular four-dimensional lattice 
and perform a Fourier transform. Now there will be four momentum variables Pi,p 2 ,P3,P4, all 
ranging from — % to it. A convenient form to explore is 

z = B(AC — cos(/?i) — cos(p2) — cos(/?3) — cos(p4)) 
+io x ( sin(pi ) + s'm(p 2 ) -sin(p 3 ) -sin(p 4 )) 
+iOy(s'm(pi) -sm(p 2 ) -sin(p 3 ) + sin(/? 4 )) 

+ia z (sm(p l )-sm(p 2 )+sm(p 3 )-sm(p 4 )). 
Here B and C are parameters whose values will be discussed later. 

For zero energy states we need z to vanish. This gives four equations corresponding to the 

coefficients of 1 and each a, vanishing. From the coefficients of the Pauli matrices imply the sines 

of all the momentum components must be equal. Picking p\ arbitrarily, each other p^ must either 

equal p\ or % — p\ . Now turning to the constant part of z, we have 

cos(pi) + cos(p 2 ) +cos(j?3) +cos(/?4) = 4C. (12) 

Since the cosine function is bounded by unity, we clearly must take C < 1 to have any solutions. 
To resolve the pi <-> 7T — pi ambiguity it is convenient to ask that cos(/?,) be positive. Imposing the 
constraint C > 1/2 ensures this. We will later discuss some interesting specific choices for C. 

So with these constraints on the constant C there are exactly two zeros of energy in the Brillouin 
zone. These occur when all components of p are equal and satisfy cos(p) = C. The two solutions 
differ in the sign of p. Picking the positive sign for convenience, it is useful to expand about the 
zero 

Pii=P + qn (13) 



with cos(j5) = C and p > 0. Defining S = sin(/5) = a/1 — C 2 , we have 

cos(du) =Ccos(^r„) -Ssin(<7„) = C-Squ +0(q 2 ) 

(14) 

sin(p M ) = Scos(q^) + Csin(q^) = S + Cq fl + 0(q 2 ). 
Inserting all this into our quaternion 

Z = BS(qi +<?2 + <?3 + <?4) 
+iCa x (#1 + q 2 - q3 - 94 ) 

(15) 

+iCo y (qi -q2~qi+ q*) 

+iCa z (qi-q 2 + q3-q4) +0(q 2 ). 
This two by two matrix is to be inserted into the analogue of Eq. © to otain a four by four matrix. 
At this point we introduce a convention for the Dirac gamma matrices 

7 = o x ® d 

y A = o y ®\ (16) 

y 5 = O z ®\ =7172/374 

The direct product notation here is defined so that 75 is diagonal with —1 in the last two places. 
With these conventions our Euclidean Dirac operator takes the form 

D = C(qi +q 2 -q3~ #4) «7l 
+C(qi -<?2-<73 + <74)z'72 

(17) 

+C(qi -qi + q-i-q^iy-h 
+BS(q l +q 2 +q3+ 94>'74 + 0(q 2 ) . 
This reproduces the desired massless Dirac equation if we identify new momenta 

k\ =C(qi+q 2 -q3-q A ) 
k 2 = C(qi-q 2 -q 3 + q4) 

(18) 

h = C(qi-q 2 + q3-q 4 ) 
k A = BS(q 1 + q 2 + q 3 + q A ) . 
Proper Lorentz invariance requires symmetry between the k's. This implies that the original lattice, 
as generated by translations using the q's, will in general be distorted from simple hyper-cubic. 

The physical angles between the original lattice directions are easily determined. For example, 
a step along the q\ axis represents a shift in the basis by (C,C,C,5S) while going along the q 2 
axis gives (C, — C, —C,BS). These vectors are at the angle 

d2 r>2 /-Q. 
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to each other. 

Note that the original axes can be made orthogonal by setting B = C/S = cot(p). With such 
a choice, gauging the theory is straightforward. Because the starting links are orthogonal, we 
can use the usual Wilson gauge action with group elements on links interacting with the simple 
plaquette action. Borici |7|] has recently suggested using 5=1 and S = C = 1/ \fl as a particularly 
simple case with orthogonal axes. 

Another especially interesting choice for the parameters B and C gives a closer analogy with 
graphene. Imagine the fermion fields on one site to actually be spread along a new bond connecting 
two atoms in the k$ direction, similar to the construction for graphene indicated in Fig. (0Q). This 
structure would be particularly symmetric if the angles between all five bonds attached to an atom 
were equal in length and distributed at a common angle to each other. 

To construct this structure it is simplest to look in momentum space continued beyond the first 
Brillouin zone. The zeros of the energy then form a lattice on which we wish to impose the above 
symmetry. The idea is to consider one zero and ask that the five closest ones are equally distant in 
the physically scaled momentum units and at angles from each other satisfying cos(<9) = —1/4. 
More precisely, ask that the physical four vector between the two zeros at = p and p^ = —p 
be of the same physical length as the vector from the first of these to the repetition of the second 
when any one of the four components of p is increased by lit. We obtain a set of five four vectors, 
which should be symmetrically distributed in four space. The conditions on lengths and angles 
then determines the parameters as 

P=\ 

C = cos(7T/5) (20) 
B = V5cot(rc/5). 

This lattice of zeros has an appealing intuitive geometric interpretation in terms of bonds along 
one direction, analogous to the horizontal bonds in Fig. CD, splitting at a site into four bonds 
going off symmetrically into four-space. These then join a repetition of this structure with new 
horizontal bonds displaced in the various directions. While in two dimensional graphene each 
carbon is coupled symmetrically to three neighbors, here each atom is directly coupled to 5 others. 
The entire lattice is then built up of hexagonal chairs with an inter-bond angle of acos( — 1/4) = 
104.4775 . . . degrees. Note that the diamond lattice in three dimensions represents an intermediate 
case, where one bond splits into three giving each atom a tetrahedral environment. 

This lattice has a 120 element discrete symmetry group under permutations of the five neigh- 
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bors. The specific action chosen here, however, selects the £4 direction as special from the way 
it appears in the site diagonal term. Because of this, the full action is not exactly invariant un- 
der the full 120 element group, but only a 24 element subgroup of rotations and inversions that 
leave invariant the positive diagonal of the initial hypercubic structure defined by the q axes. This 
tetrahedral symmetry, which is that of the three dimensional diamond lattice, applies for all valid 
values of the parameters B and C and includes 12 parity odd operations. Considering the time axis 
to run along this positive diagonal, there is in addition a symmetry under reflections about a space- 
like surface perpendicular to this time axis accompanied by a rotation to bring the bonds back into 
alignment. These symmetries will eliminate many possible lattice artifacts and constrain operator 
mixing in renormalization studies |8|]. The main remaining issue is that the special treatment of 
time still allows a renormalization of the speed of light. This can be compensated by adding to the 
gauge action terms involving six link chairs orthogonal to the time direction. While such terms are 
generally necessary, for the hyper-cubic values B = C/S they should go to zero in the continuum 
limit. Indeed, these are all issues of order the lattice spacing squared, as discussed in Q9[| . 

From the standpoint of computational efficiency, it does not appear to matter much what val- 
ues of the parameters B and C are chosen, within the constraint 1/2 < C < 1. Near the limits 
of this range one should expect lattice artifacts to increase. The hyper-cubic values B = C/S are 
presumably closest to conventional lattice gauge ideas and simplest to gauge, with only a tradi- 
tional plaquette action required. On the other hand, the values in Eq. (|20l) may have smaller lattice 
artifacts due to the high symmetry. 

Returning from reciprocal to position space, the fermionic action involves several terms. First 
from the collapsed £4 bond there is a site diagonal term AiBCxyy^xir. Then for a forward step 
in the various directions we pick up a factor of the hopping parameter K multiplied by different 
combinations of gamma matrices, as listed here: 

for a hop in direction 1 : +ji + 72 + 73 — iBf4 
for a hop in direction 2: +Y\ — 7? — 73 — iByk 

(21) 

for a hop in direction 3: —J\ — 72 + 73 — Z.B74 

for a hop in direction 4: — 71 + 72 — 73 — iBj4 
Keeping the operator D anti-Hermitean, the reverse hops involve minus the conjugate of these 
factors. Note the factor of i in front of 74 which is absent for the 7i_3 terms. This twisting of the 
phase gives rise to the required factors of sin or cos in Eq. (flTT) . This action is only marginally 
more complicated than that of naive fermions; so, it should be easy to insert into simulations. 
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Chiral symmetry is manifested in the exact anti-commutation of 75 with D. This is actually a 
flavored chiral symmetry since the expansion about the negative solution for p flips the sign of the 
gamma matrix associated with ]cq. Note that as with naive fermions, D is purely anti-Hermitean. 
The chiral symmetry can easily be broken with the addition of a term proportional to 75 to H that 
splits the degeneracy of the type a and type b sites. This gives each physical fermion a common 
mass. Because the chiral symmetry remains exact on gauging, an additive renormalization for the 
physical fermion mass is forbidden, unlike with Wilson fermions. The further addition of a Wilson 
type mass term would enable splitting the degeneracy of the two species. Note that if we wish to 
extend this formalism to more flavors, the no-go theorem restricts us to an even number of species. 
The best we can do for three flavors is to start with four and, using a chiral symmetry breaking 
operator, make one of the flavors heavier. Such a term will unfortunately introduce the possibility 
of an additive mass renormalization for the third flavor. 



It is perhaps interesting to compare this formalism with staggered fermions ifioi fll. 



jj . The 



latter also have an exact chiral symmetry and the appearance of several degenerate flavors (some- 
times called tastes to distinguish them from independent lattice fields). In both formalisms the 
chiral symmetry is flavored, with the different species rotating differently under a chiral rotation. 
Also, in both cases mixing between the flavors will result in the zero modes associated with topo- 
logically non-trivial gauge fields to no longer be exact at finite lattice spacing IU3U 14I1. However 
the usual staggered approach has the doublers appearing in multiples of four, and thus is further 
from the situation desired for the physical case with only the light up and down quarks. Staggered 
fermions do have only one component of the field per site, yielding potentially faster simulations 
at the expense of an intricate non-local formalism in constructing composite fields. 
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